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mean  of  the  power.  An  expression  is  obtained  for  the  probability  distribution 
bfthe  amplitude  and  power  of  the  error  field  assuming  the  error-free  field  is 
real  and  there  are  sufficient  elements  in  the  array  for  the  Central  Limit 
Theorem  to  be  applied.  For  small  phase  errors,  and  points  not  halfway 
between  grating  lobes,  the  expression  reduces  to  the  Rice-Nakagami  distribu¬ 
tion.  It  further  reduces  to  the  Rayleigh  distribution  at  points  where  the  error- 
free  pattern  has  a  null.  Expressions  for  the  pattern  power  mean  and  variance 
for  an  array  of  arbitrary  size  obtained  directly  are  consistent  with  those 
obtained  by  assuming  a  large  number  of  array  elements  and  using  the  Central 
Limit  Theorem  to  derive  the  probability  distribution  of  the  pattern  power. 
-..—ik.The  theoretical  results  are  compared  with  Monte  Carlo  computer  simula¬ 
tions.  The  Kolmogorov-Smirnov  test  usually  shows  close  agreement  between 
the  simulation  and  the  theoretical  results.  However,  the  test  also  shows  that 
the  depths  of  multiple  nulls  are  not  statistically  independent  when  there  are 
multiple  nulls  in  the  error-free  pattern;  the  probability  distribution  of  the  depth 
of  the  least  deep  null  is  closer  to  the  distribution  predicted  for  a  single  null 
than  it  is  to  the  distribution  for  multiple  nulls. 
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Statistical  Analysis  of  the  Effect  of  Phase 
Quantization  on  Array  Antenna  Sidelobes 


1.  INTRODUCTION 

There  is  currently  considerable  interest  in  using  phased  arrays  with  digital 
phase  shifters  for  adaptive  nulling.  However,  the  quantized  output  of  a  digital  phase 
shifter  differs  from  the  ideal  value.  This  difference  between  the  actual  and  ideal 
phases  is  reflected  in  a  corresponding  discrepancy  between  the  actual  pattern  and 
the  ideal  or  error-free  pattern.  Of  special  concern  for  adaptive  nulling  applica¬ 
tions  is  the  fact  that  phase  quantization  errors  can  cause  shallower  nulls,  reducing 
the  ability  of  the  array  to  reject  interference  from  particular  directions.  There¬ 
fore,  understanding  the  influence  of  phase  quantization  errors  on  array  antenna 
patterns  is  important. 

A  standard  technique  used  to  eliminate  the  spurious  quantization  lobes,  similar 
to  grating  lobes,  that  result  from  a  periodic  phase  quantization  error  distribution 
across  the  elements  of  an  array,  is  to  randomize  the  phase  shifters  by  inserting 
a  constant  phase  shift  which  differs  from  element  to  element.  1  The  phase  quantiza¬ 
tion  errors  can  then  be  regarded  as  independent  of  each  other.  Any  general  study 
of  the  effect  of  such  phase  quantization  errors  on  array  patterns  is  necessarily 
statistical  in  nature,  focusing  on  statistical  parameters  of  the  phase-error  induced 
pattern  perturbations  such  as  the  mean  and  variance  of  the  power  at  a  specified 

(Received  for  publication  7  July  1982) 

1.  Skolnik,  M.  I.  (1980)  Introduction  to  Radar  Systems,  McGraw-Hill,  New  York,  p,  321. 


5 


pattern  location.  Since  the  magnitude  of  the  phase  quantization  errors  has  a  maxi- 
N 

mum  value  of  til  2  Dlt,  where  N^t  is  the  number  of  bits  in  the  phase  shifters,  it  is 
natural  to  base  a  statistical  description  of  the  pattern  perturbations  on  the  assump¬ 
tion  that  any  particular  set  of  element  weight  phase  errors  is  obtained  by  drawing 
each  phase  error  independently  (with  replacement)  from  an  underlying  population  of 

Nbit  N 

phase  errors  uniformly  distributed  in  the  interval  [-7t/2  ,  til  2  ]. 

This  report  describes  the  statistical  properties  of  the  sidelobes  of  linear  arrays 

whose  element  phases  are  independently  and  uniformly  distributed.  The  array 

elements  are  assumed  to  be  isotropic  and  equispaced  and  to  have  a  deterministic 

but  arbitrary  amplitude  taper.  A  number  of  investigations  have  been  devoted  to 

2 

various  related  aspects  of  this  subject,  and  in  considerable  part  this  report  is 
designed  to  bring  together,  for  convenient  reference  and  in  self-contained  form, 
results  that  are  well  established  in  the  literature. 

We  begin  by  investigating  tne  mean  and  variance  of  the  power  at  an  arbitrary 
pattern  location.  No  conditions  are  imposed  regarding  the  number  of  array  ele¬ 
ments,  The  resulting  expression  obtained  for  the  variance  of  the  power  is  believed 
to  be  new.  Approximations  valid  for  small  phase  errors  (N^  >4)  are  obtained  from 
the  general  expressions  for  the  mean  and  variance.  We  then  derive  an  expression 
for  the  probability  distribution  of  the  amplitude  and  power  of  the  perturbed  field 
under  the  assumptions  that  the  error-free  field  is  real  and  that  the  number  of 
elements  in  the  array  is  sufficiently  large  for  the  Central  Limit  Theorem  (CLT)  to 
be  applied.  We  specialize  this  distribution,  first  to  the  case  of  small  phase  errors 
and  ordinary  sidelobe  locations  (locations  not  halfway  between  grating  lobes),  and 
then  to  the  important  case  of  sidelobe  locations  for  which  the  error-free  pattern 
has  a  null.  We  show  that  under  the  conditions  for  which  the  CLT  is  valid,  the  exact 
expression  for  the  variance  of  the  power  obtained  with  no  restriction  on  the  number 
of  elements  reduces  to  the  expression  obtained  by  means  of  the  CLT  and  the  prob¬ 
ability  distribution  of  the  power. 

The  statistical  distribution  of  power  at  the  least  deep  null  among  several  loca¬ 
tions  for  which  the  error-free  pattern  has  a  null  is  then  obtained  by  assuming 
independence  of  the  null  depths  at  two  or  more  locations  and  multiplying  together 
the  probability  distributions  for  the  individual  nulls.  The  theoretical  analysis  is 

2.  See  for  example,  Ruze,  J.  (1952)  The  effect  of  aperture  errors  on  the  antenna 

radiation  pattern,  Nuovo  Cimento  Suppi.  ,  ...Si  No.  3):364-380;  Allen,  J.  L.  (1961) 
Some  extensions  of  the  theory  of  random  error  effects  on  array  patterns. 

Chap.  Ill,  Phased  Array  Radar  Studies,  1  July  1960  to  1  July  1961,  Al.  1.  T. 
Lincoln  Lab  Tech.  Rep.  236;  Skolnik,  M.  I.  (1980)  Nonuniform  arrays. 

Chap.  6  of  Antenna  Theory,  Part  1,  Collin,  R.  E.  ,  and  Zucker,  F.  J.  (Eds.) 
McGraw-Hill,  M.  Y.  ,  pp.  318-321;  Steinberg,  B.  D.  (1976)  Principles  of 
Aperture  and  Array  System  Design,  Chaps.  8,  9,  and  13,  John  Wiley,  N.  S'., 
Shifrin,  Y.S.  ( 197  1)  Statistical  Antenna  Theory,  Chap.  7,  Golem  Press, 
Boulder,  Colorado. 


concluded  by  deriving  the  distribution  of  power  at  locations  halfway  between  grating 
lobes.  Again  we  demonstrate  consistency  between  the  exact  expression  for  the 
variance  of  the  power  and  the  expression  derived  from  the  CLT  under  the  assump¬ 
tion  of  large  arrays.  Following  the  theoretical  analysis,  we  describe  the  results 
of  computer  simulations  performed  to  compare  with  theoretical  results. 


2.  ANALYSIS 

In  this  section  we  derive  expressions  for  the  mean  and  the  variance  of  the 
power,  and  the  probability  distributions  of  the  amplitude  and  power,  at  a  given 
direction  of  the  pattern  of  a  linear  array  of  equispaced  isotropic  radiators,  given 
that  the  phases  of  the  array  elements  are  subject  to  random  errors.  To  investigate 
the  effect  of  phase  shifter  quantization  on  the  array  pattern,  we  take  the  mean, 
variance,  and  probability  distribution  over  an  ensemble  of  arrays;  in  each  member 
of  the  ensemble,  we  assume  the  phases  of  the  elements  to  have  mutually  independent 
random  deviations  6  from  their  ideal  values,  uniformly  distributed  in  the  interval 
r  N  j  .  N  b  •  ."I 

J^-jt/2  n  ,  t;  I  2  1  J  where  N^-t  is  the  number  of  bits  in  the  phase  shifters. 

Let  d  be  the  spacing  between  the  array  elements  and  assume  the  phase  refer¬ 
ence  center  to  be  the  center  of  the  array.  Let  w^  be  the  ideal  complex  weight 
(that  is,  the  weight  in  the  absence  of  phase  shifter  errors)  of  the  nth  element,  and 
let  6^  be  the  corresponding  random  phase  error.  Then  the  error-free  field  pattern  is 


.  Y*  jdn' 

0  =  w  e 


md  the  pattern  in  the  presence  of  phase  errors  is 

Z1  j  6  j  d  u 

w  e  ne  n 

n 


d  =  -  (n  -  1),  n=  1,  2, _ N 

n  2 


u  =  kd  sin  0 
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k  =  2  tt/X 


and  8  the  angle  measured  from  broadside  to  the  array. 


2.  L  Mean  of  the  Power 


We  first  obtain  an  expression  for  the  ensemble  average  of  the  power  at  a  given 
direction  of  the  pattern.  The  pattern  power  of  a  member  array  is  given  by 


I  F(u)  | 


=  F(u)F(u)  =  w  w  e  n 

“  Tr*'  n  m 


1(6-6  )  j(d  -d  )u 

J  n  r-n  J  n  rr\ 


m 


(2) 


n  m 


so  that  we  seek  the  ensemble  average,  |F(u)|  .  In  calculating  |F(u)|  we  make 


use  of  the  theorem  that  the  expectation  of  a  product  of  mutually  independent  random 


variables  is  equal  to  the  product  of  their  expectations.  Hence,  in  the  double 


summation  in  Eq.  (2)  we  separate  the  terms  for  which  n^m  and  hence  6n  and  6^ 


are  independent)  from  those  for  which  n=  m  (and  hence  6=6).  For  n=  m,  the 

n  m 


contribution  to  the  ensemble  average  is  simply 


while  for  n^  m  we  have 


y-  y  w  *  "j(V6^  pj<Vdm)u 
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* — '  A —  n  m 

n  m 

n^  m 


J5 


?  £  % 
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n  m 
n4  ni 


w  e 
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/here  we  have  made  use  of  the  fact  that 


j(8  -  5  )  ITT 

J  n  m  J  n 
s  =  e 


=  |ejd 


N. 


Nl 


r  bit~l 

with  6  a  random  variable  uniformly  distributed  in  I-tt/2  ,  tt 1 2  |  since  the 


tl  . 

sir 


phase  errors  are  assumed  to  be  identically  randomly  distributed  for  all  elements 
of  the  array.  The  ensemble  average  of  the  power  is  then 


3.  Feller,  W.  (1968)  An  Introduction  to  Probability  Theory  and  Its  Applications, 
Vol.  1.  3rd  Ed.  .  John  Wilev,  N.  Y.  ,  pp.  222  -230. 


i’<u) 


|wn|2  + 


7F|2  V'  V'  . *  .j(Vdm)u 


S  E  % 


n  m 
n4  m 


w  e 
m 


z  z 

n  m 


w  w  e 
n  m 


i(d  -  d  )u 
J  n  m 


+  (1 


^l2.  z 


But 


z  z 

n  ni 


w  w  e 
n  m 
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n  m  =  Ko(u)Fo'(u)  =  |Fo(u)| 


so  that 


F(u) 


"J7 
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The  calculation  of  e-1  is  straightforward.  For  phase  errors  uniformly  dis¬ 


tributed  in  the  interval  [  -A,  A|,  A  =  7t/2 


bit 


=  23  f  °j6  d6 

-A 


sin  A 

s~  • 


(3) 


Hence 

pr=(^)2 

and 


|  F(u) 

The  ensemble  average  power  pattern  is  thus  the  superposition  of  two  terms,  the 
first  of  which  is  the  error-free  power  pattern  multiplied  by  a  factor  whose  magni¬ 
tude  is  less  than  unity  for  A  >0,  and  the  second  is  independent  of  the  direction. 


(sin  A  1 


F  (u) 
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H^)1  z 


(4) 
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sin  A 


For  A  =  0, 


shifters  with  nuinv  bits  (N,  ..  >  4), 

bit 


1  and  I  F(u)|  “  reduces  to  |  F  (u)|  In  applications  to  phase 


()■>  ■> 

A“ 

=  '~T 


and 


F(u> 


=  ('--) 


(u) 


^  E  iwj2 
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At  a  null  of  the  error-free  pattern 
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!  F(u)  | 
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A" 
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/  sin  A  \ ‘ 

\  A  / 


E 


E  w  ,  !•  <«)  -  0 

n  o 
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N. >  4  • 
bit 


(6a) 


(6b) 


Note  that  for  each  extra  bit  of  the  phase  shifters,  A“  is  multiplied  by  a  factor  of  1/4 
and  hence  the  average  power  at  a  null  of  the  error-free  pattern  decreases  by  6  dB. 
The  significance  of  the  filling-in  of  the  nulls  of  the  error-free  pattern  due  to  phase 
errors  is  usefully  assessed  in  relation  to  the  peak  value  of  the  error-free  pattern. 
For  error  free  weights  such  that  the  phase  varies  linearly  across  the  array,  the 


peak  power 
obtain 


Dividing  Eq.  (6a)  by  this  quantity  to  normalize  it,  we 


F  ( u )  1 


(?" 
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2.2  Variance  of  the  Power 


’'•I 


Next  we  obtain  an  expression  for  the  variance  of  the  power 


Since  we  have  already  found  |  !•*( n ) I  .  we  must  determine 


[ i  Kuo  i  “  L  E  E  E w  w m  « 

L  n  m  i)  (i 


-  (1  4  0  -01  |((1  -  (I  4  <1  -  (I 

*  •  n  m  p  (j  n  in  p 

we  '  1  e  1 
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As  above  in  obtaining  an  expression  for  the  average  power,  we  proceed  by  grouping 
the  terms  of  the  quadruple  summation  according  to  whether  or  not  equalities  hold 
among  the  indices.  The  following  cases  must  be  distinguished: 

(1)  All  four  indices  are  equal 
n=m-p-q 

(2)  three  indices  equal 
(2a)  n-m-pVq 

(2b)  n-m-q^p 
(2c)  n=p-q^m 
(2d)  m-p-q#n 

(.1)  two  pairs  of  indices  equal 
(da)  n-  m,  p-q,  n*p 
(dh)  n=p,  m-q,  iWm 
(dc)  n~ q,  ni-p,  iW in 

(4)  two  indices  equal 

(4a)  n-  ill,  tWp,  tWq,  p^q 
(lb)  n-p.  iWm,  n^q,  mV <) 

(4c)  n-q,  ttfm,  n*p,  m^p 
( 4 cl )  ni-p,  lilVn,  m^q,  n?q 
(4e)  m-q,  m^n,  m*p,  n^p 
( 4 f )  p-q,  p^n,  p-lm.  Him 

(a )  no  indices  etpial 

We  treat  each  of  these  cases  in  turn. 

(  1 )  n-  m-  p-q 


The  contribution  of  the  corresponding  terms  in  the  quadruple  sum  of  Eq.  (8) 
is  simply 


(2a)  n=m=p^q. 

The  contribution  of  the  appropriate  terms  in  the  quadruple  sum  is 
«— *  '■>  ■  j(6  -  (j  )  j(d  -d  )ti 
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where  Gq(u)  is  the  pattern  corresponding  to  the  weights  |  w^l  w^;  that  is, 
N  .  . 

Eq  1  d  U 

|  w  I  w  e  11 
n=  1  n  11 


(2b)  n=ni=q^p 

Here  the  contribution  of  the  quadruple  sum  is 
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(u)  1-  (u) 


(2c)  n=p=q^m 

By  symmetry,  the  contribution  of  these  terms  is  the  same  as  that  for  (2a). 
(2d)  m=p=q^n 

By  symmetry,  the  contribution  of  these  terms  is  the  same  as  that  for  (2b). 
(3a)  n=m,  p=q,  n^p 

The  contribution  of  the  appropriate  terms  in  Eq.  (8)  is 
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(3b)  n=p,  m=q,  n=^m 

The  contribution  of  these  terms  is 
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where 


\ 

Ho(u)  =  Yj 

n=  1 
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n 


1  2  d  u 
J  n 
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(10) 


Note  that  H  (u)  is  the  pattern  of  an  array  whose  weights  are  the  squares  of  the  given 
array  and  whose  elements  are  spaced  twice  as  far  apart. 
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(3c)  mq,  m=  p.  n^m 

Bv  symmetry,  the  contribution  of  these  terms  is  the  same  as  that  for  (3a). 
(4a)  n  m,  n/p,  n/q,  p^q 

The  corresponding  terms  in  Eq.  (8)  contribute 
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The  eontribution  of  these  terms  is 
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where  Gq(u)  is  given  by  Eq.  (9)  and  Hq(u)  by  Eq.  (10). 


(4c)  n=q,  n^m,  nyp,  m^p 

By  symmetry,  the  contribution  of  these  terms  equals  that  of  (4a). 

(4d)  m=p,  m^n,  m^q,  n^q 

By  symmetry,  the  contribution  of  these  terms  equals  that  of  (4a). 

(4e)  m=q,  m^n,  m^p,  n^p 

By  symmetry,  the  contribution  of  these  terms  equals  the  complex  conjugate 
of  that  of  (4b). 

(4f)  p=q,  p^n,  p^m,  n^m 

By  symmetry,  the  contribution  of  these  terms  equals  that  of  (4a).  Hence  of 
the  six  cases  for  which  exactly  two  indices  are  equal,  four  are  equivalent  to  (4a) 
and  two  are  complex  conjugates. 

(5)  no  indices  equal 

These  terms  contribute 
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To  obtain  the  first  equality  we  have  simply  pooled  the  results  of  the  preceeding 

cases.  _ 

The  average  of  the  squared  power,  |  F(u)|  ,  is  equal  to  the  sum  of  the  con¬ 

tributions  of  the  separate  cases  we  have  analyzed.  Equation  (7)  shows  that  to 
calculate  the  variance  of  the  power  we  must  subtract  the  square  of  the  average 
power  from  the  average  of  the  squared  power.  From  Eq.  (4)  the  square  of  the 
average  power  is  given  by 
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Subtracting  Eq.  (11)  from  |  F(u)|  4  and  combining  terms,  the  variance  of  the  power 
is  found  to  be 
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We  have  shown  above,  in  Eq.  (3),  that  for  phase  errors  uniformly  distributed  in  the 
interval  | -A,  A),  A  =  n/ 2^bi\  ^  .  Similarly  J77  =  ***** 


and  so 
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where  G^iu)  and  iMu)  are  given  by  Eqs.  (9)  and  (10)  respectively.  As  above  for  the 
average  power  pattern,  the  variance  of  the  power  is  the  superposition  of  terms 
that  are  direction  dependent  and  terms  that  are  direction  independent  (that  is, 
omnidirectional).  For  A  =  0  (that  is,  no  phase  errors),  ^  and  are  both 

unity  and  Eq.  (12)  reduces  to  zero  as  required.  In  applications  to  phase  shifters 
with  many  bits  we  can  approximate  the  combinations  of  and  that 

appear  in  Eq.  (12)  by  their  small  angle  forms  thus  obtaining  after  some  manipula¬ 


tor  ^  |  F(u)|  )  = 
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It  is  of  interest  to  estimate  the  relative  contribution  of  the  three  terms  within  the 

square  brackets  on  the  left  hand  side  ofEq.  (14).  For  phases  of  the  error -free  weights 

-j  d  u . 

varving  linearly  over  the  array  we  can  write  =  a^e  n  b  and  so  from  Eq.  (10) 
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Hence  H  (u)  has  a  maximum  magnitude  of  a2  at  those  directions  for  which 
o  n 

u-  u  t  m  n  or  equivalently 

k  cl  sin  ft  -  k  d  sin  0  i  mi 


sin  0  =  sin  ft  ±  m  __ 
s  771 


If.  for  exant pie,  d./.'=  1  and  ft  =  0  then  in  the  visible  region  Hq(u)  =  a2  when 
sin  fi  =  0,  :  1/2  or  i)  0°,  i  30°.  As  u  moves  away  from  these  points,  H  (u)  fall; 
off  rapidlv,  especially  for  arrays  with  a  strong  taper  of  the  amplitudes  of  the 


error-free  weights,  and  so  ill  (u) 
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As  far  as  the  relative  contribution  of  the  terms  ^  and  §  I  w„ 

concerned,  applying  Cauchy's  inequality  (Ref,  4)  we  obtain 

E  u.l*  ■  E  [(kifl  *  ( E  i./)2 

n  nLN  'A  n  ' 

and  applying  Chebyshev's  inequality  (Ref.  5)  we  obtain 

4 1-.'4-  k  [M)2m?ki2)2 

so  that 


Note  that  the  left  "s"  in  Hq.  (15)  becomes  equality  for  the  case  of  a  uniform  ampli¬ 
tude  array,  and  that  the  right  "s"  becomes  equality  for  the  case  of  an  array  with 
all  amplitudes  except  for  one  equal  to  zero.  Furthermore  it  can  be  shown  (see 
Appendix  A)  that  for  any  amplitude  taper  describable  by  a  polynominal  in  n,  the 


4.  Cauchv's  Inequality; 

(equality  for  a.  -  c  h 

K  K 

e  constant). 

See  Abramowitz,  M.  (1P72)  Hlemcntarv  analytical  methods.  Chap.  5  of  Handbook  of 
Mathematical  Functions.  M.  Abramowitz  and  I.  Stegun,  (Kds.  ),  I )over,  \.  V. 

5.  Chebvshev's  Inequality; 
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See  Abramowitz,  M.  . 
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ratio  of  7  ,  |  w  1  to  ZJ  w  |  J  is  proportional  to  1/N  as  N  becomes  large.  Hence 

6  v**^  i  i  4 

tor  large  arrays  and  any  usual  amplitude  taper,  ■=■  /  ^  I  w  I  is  of  the  order  of  1/N 

(  2  3  n 
I  w  |  j  and  can  be  neglected. 

It  is  often  useful  to  express  the  variance  of  a  quantity  relative  to  the  square  of 


compared  with 

I2"  I  w^|  j  and  can  be  neglected. 

It  is  often  useful  to  exfpress  the  variance  of  a  quantity  relative  to  the  square  of 
the  mean  of  the  quantity,  thus  giving  a  measure  of  variance  divorced  from  the  actual 
units  in  which  the  quantity  is  measured.  [The  ratio  of  the  standard  deviation  (that 
is,  the  square  root  of  the  variance)  to  the  mean  is  known  as  the  coefficient  of  varia¬ 
tion.  ]  If  Eq.  (14)  is  divided  by  the  square  of  Eq.  (6b)  we  obtain 


Note  that,  to  within  the  restriction  that  A  is  small  enough  for  (sin  A)/A  and 
(sin  2  A)/ (2  A)  to  be  well  approximated  by  their  small  angle  forms,  the  ratio  of  the 
variance  of  the  power  to  the  square  of  the  mean  power  at  a  direction  for  which  the 
error-free  pattern  has  a  null  is  independent  of  the  magnitude  of  the  phase  errors. 

In  accordance  with  the  discussion  of  the  previous  paragraph,  if  u  is  not  in  the 
vicinity  of  those  points  for  which  Hq(u)  has  its  maxima,  and  if  the  number  of 
elements  in  the  array  is  large,  the  second  term  on  the  right  hand  side  of  Eq.  (16) 
can  be  neglected  compared  with  one  and  so  the  variance  of  the  power  is  approxi¬ 
mately  equal  to  the  square  of  the  mean  of  the  power.  Equivalently,  the  standard 
deviation  of  the  power  is  approximately  equal  to  the  mean  of  the  power. 

2.3  Probability  Distribution  of  the  Power 

Having  derived  expressions  for  the  mean  and  variance  of  the  power,  we  now 
obtain  the  probability  distribution  of  the  amplitude  of  the  field  in  a  given  direction 
of  the  pattern.  By  a  simple  transformation,  we  can  then  also  obtain  the  probability 
distribution  of  the  power  in  a  given  direction.  Our  procedure  is  essentially  that  of 
Allen  and  Beckmann.®  We  first  resolve  the  field  into  its  real  and  imaginary 
components,  each  of  which  is  given  by  a  sum  of  N  terms  (N  =  number  of  elements 
in  the  array)  involving  trigonometric  functions  of  the  random  phase  errors.  Since 
the  phase  errors  are  assumed  to  be  mutually  independent,  the  terms  of  the  sums 
are  independent  random  variables.  We  assume  that  N  is  sufficiently  large  (and  that 


6.  Beckmann,  P.  (1967)  Probability  in  Communication  Engineering.  Harcourt. 
Brace  and  World,  Inc.  ,~New  York.  op.  S9-63. - S - 


no  one  term  dominates  the  sum)  for  the  Central  Limit  Theorem  to  apply,  and  so 
infer  that  the  real  and  imaginary  parts  of  the  field  are  Gaussian  distributed.  It  is 
straightforward  to  calculate  the  mean,  variance,  and  covariance  of  the  real  and 
imaginary  field  components.  We  then  transform  from  cartesian  to  polar  coor¬ 
dinates  and  obtain  the  probability  distribution  of  the  amplitude  of  the  field  by  inte¬ 
grating  the  joint  probability  distribution  of  the  amplitude  and  phase  of  the  field  with 
respect  to  the  phase. 


Referring  to  Eq.  (1),  let  the  error-free  weights  be  represented  by  a^  e 
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Then 
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and 
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Since  our  interest  is  in  large  arrays,  and  since  for  such  arrays  no  one  element  can 
be  singled  out  as  dominant,  we  can  apply  the  Central  Limit  Theorem  to  assert  that 
the  real  and  imaginary  parts  of  the  field  are  random  variables  with  a  joint  bivariate 
Gaussian  probability  distribution.  The  statistics  of  this  distribution  are  completely 
determined  by  the  means,  variances,  and  covariance,  so  that  our  task  is  to  calcu¬ 
late  these  quantities. 


where  we  nave  made  use  of  the  trigonometric  identities 


cos2  8  =  ( 1  +  cos  2  8  ) 


(22) 


sin~  8  =  ( 1  -  cos  2  8  ) 


Hence 


_2  1  f  /  sin  A  \  “  1  2  1  [7sin  A  \  2  sin  2A"I  2  , 

ax/  2  [  1  -  1^-A-j  Jan‘  2  [(-S-j  TK  J  anCOS<2^n 

_2  _  1  P  /  sin  A  \  2  1  2  if/  sin  A  \  2  sin  2  A  1  2 

Vn'  2  [  -  l  A  j  J  an+  2  [[  A  )  2  A  J  dn 


cos  (2 i^n> 


and 


4 


4- 


i[ 


l  - 


sin  A 


sin  2A 

2  A'  " 


sin  2  A 

2  A 


]£4 

-*  n 

]2-J 


cos  (2^n> 

(23) 
cos  (2^n> 

(24) 


The  covariance  of  the  real  and  imaginary  parts  of  the  field  is  defined  by 


axv  =  XY-XT. 


Now,  from  Eqs.  (17b)  and  (18b), 


XV  = 


E  E 

n  m 
nA  m 


a  (cos  ^  cos  6  -  sin  sin  6  )  (sin  li/  cos  6  +  cos  \1 >  sin' 

n  m  n  n  n  n  ^  m  m 


+  1  Ea„  sin[2(i//n  4  5 n )  1  (25) 

n 

where,  as  above,  we  have  separated  the  double  summation  into  the  terms  for  which 
n4  m  and  hence  6n  and  6  are  independent,  and  the  terms  for  which  n=  m.  The 
double  summation  in  Eq.  (25)  reduces  to 


I  — r — I  i__#  a  a  cos  \L  sin  t// 

l  A  /  n  ni  ’n  r 

V  /  n  m 


because  of  Eq.  (19),  and 


1  v  2--—i 
?  Z-  an  am  | 


5  +  6  >  =  Tj  £  a2  sin(2t//  ) 

n  n  2  n  1 
n 


n  +  cos  (2  \p  )  sin  (2  6  )' 
n  n  n 


1  sin  2A  V  2  . 

o  *  yr  s  si 


a  sin  (2  li/  )  . 
n  n 


Hence 


=  /  — \  I" £  a  a  c°s  p  sin  p  ~  a^ cos  sm  p  1 

\  A  /  L‘ n  m  n  m  n  m  n  "  n  ”J 


1  sin  2  A 

2  2  A 


£  an  sin  (2  *„> 


S  [(^)2  -T^]  E  •><**„> 
*  (t)'  [>.<••] 


1  [/sin  a\  2  sin  2  A]  V*  2  •  to  i  \ 

-  -  2  [[-*-)  '  ~2A~J  £  an  sin  2  ^n  * 


We  have  thus  obtained  expressions  for  all  the  quantities —means,  variances, 
and  covariance —needed  to  define  the  joint  probability  distribution  of  the  real  and 
imaginary  parts  of  the  field.  Although  it  is  possible  to  proceed  to  obtain  the  dis¬ 
tribution  of  the  amplitude  of  the  error  field  for  the  general  case  of  a  complex 

g 

error-free  field,  for  our  purposes  here  it  is  sufficient  to  restrict  attention  to  the 
case  of  a  real  error-free  field.  If  Im[  Pq(u ) J  -  0,  we  see  from  Kq.  (21)thatY=0. 
Furthermore,  if  the  imaginary  part  of  the  field  is  zero,  then  the  amplitudes  of  the 
array  element  weights  must  be  even-svmmetric  and  the  phases  of  the  weights  odd- 


symmetric  with  respect  to  the  center  of  the  array.  Hence  the  covariance  of  the 
real  and  imaginary  parts  of  the  field  vanishes.  We  then  have 


X  =  F  (u) 

A  o 


(26) 


T-  0 
_2 


(27) 


'X 


1 

=  2 

.  /  sin  a'' 

L 

i2]?^ 

i 

■  i 

[(^)2 

sin  2  A 

- J 

E 

n 

1 

r>  -(^ 

\2lv  a2 

i 

/ sin  A  \2 

sin  2  A  | 

E 

=  2 

)  J?  a- 

2 

■  »  / 

n 

cos  (2  U/  ) 
n 


(28) 


cos  (21/0 


(29) 


aXY  "  ° 


and  the  joint  probability  distribution  of  the  real  and  imaginary  parts  of  the  field  is 
given  by 


P<X.Y>  =  Wir  -P  -Tj-l  . 

7 r  ~  ”  L  2<7„  2(7,/  J 


(30; 


X  Y 


X 


We  now  transform  to  polar  coordinates  letting 
X  =  I K  |  cos  0 
Y  =  |  f|  sin  9 
dXdY  =  |f|  d|p|  df) 

so  that  |  F |  is  the  amplitude  and  9  the  phase  of  the  field.  Note  that  since 
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1 2  2  2 
F|  =  X  +  Y  ,  the  mean  power  is  given  by 


— T2  ~~2  ~~7 

F  ■  X  +  Y 


_  +  n  ^  +  n  2 

X  Y 


7.  Oppenheim,  A.Y.,  and  Schafer,  R.  W.  (1975)  Digital  Signal  Processing, 
Prentice -Hall,  New  .lersev,  pp.  24-26. 


and  hence  from  Eqs.  (20),  (21),  (23).  and  (24).  we  have  in  general  (that  is.  with 

no  restriction  to  real  error-free  fields) 

1^=  (^)2  |Fo«»»i2*  [>-(^)2]E  4 

in  agreement  with  Eq.  (4)  which  gives  the  average  power  obtained  more  directly 
above.  Substituting  Eq.  (31)  into  Eq.  (30). 


2n  axoy 


( I  F  |  cos  0  -  X)2 


(If  I  sin  0)' 


=  I  F  I 

-  2 n  axay 


Using  the  trigonometric  identities  of  Eq.  (22)  we  obtain 


.  -  I  F  I  -D  -P  cos  (2  0  )  +  Q  cos  0 
Zlt  CTX°\ 


whe  re 


(32a) 


(32b) 


(32c) 
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The  probability  distribution  of  the  field  amplitude  is  then  found  by  integrating 
p(  |  F 1 .  9  )  with  respect  to  0  . 


p(  I  1 ' ! )  =  (  ( i 1 ' 1 .  9  >  d  6 

J0 

=  2j'  p(  hi,  9  )  cl  0 


"  a\a\ 


j  ^-P  cos  (2  9  )  +  Q  cos  0  y 

Jo 


IJ roceeding  us  in  Kefe fence  B  we  obtain  an  expression  for  the  integral  1  in  terms  of 
a  series  of  modified  Bessel  functions,  lirst,  using  the  identities  in  deference  8 


^  E 

k=  0 


tk  lk(/.)  cos  tkfli,  ik  =  { 


1.  k-  0 


2,  k  2  0 


1,  (-/.)=  ( -  1 )  [  (/.) 
k  k 


we  call  w  riti 


-BeOS  (20) 


*  E 


(-  1 )  '  t .  I  (P)  cos  (2  k  (?) 
k  k 


ii 

i  -■  ^  (-l)k  t.  I.  (P)  (  C,Q  cos 9  cos(2ke) 

n  k  k  Jn 


Then,  using  the  relation  in  deference  8 


inU)  = 


7.  COS  9  ,  „  ,  ,  „ 

e  cos  (n  6  )  d  9 


8.  Olver,  K.  W.  J.  (1972)  Bessel  functions  of  integer  order.  Chap.  9  in  Handbook 

of  Mathematical  Functions,  M.  Abramowitz  and  I.  Stegun,  Eds.  ,  Dover,  N.  V 


we  have 


44 

X 


eQcOS  6  cos(2k0)  dS  =  ^  I9k(Q) 


so  that 


I  =  n  Z  (-l)k  t.  1.  (F)  1...  IQ) 
k=0  k  k 


and  the  desired  probability  distribution  of  the  field  amplitude  is 


p(  1 1'  I )  = 


F  -D 


axav 


X  (-1) 


k  =  0 


‘k  VP) 


(33) 


with  D,  P,  and  Q  given  bv  Kqs.  (32a-c).  The  distribution  can  be  normalized  by  letting 


r2  X2 

H  -  - Tf - 2 

+cjv 


(34a) 


If!  = 


Vax2+v 


(34b) 


and 


K  = 


ffY 


(34c) 


2 

so  that  B  is  the  ratio  of  the  constant  component  of  the  power  to  the  average  power 
in  the  fluctuating  component  of  the  power,  and  K  is  a  measure  of  the  asymmetry  in 
the  average  proportion  of  the  fluctuating  component  of  the  power  associated  with 
the  real  and  imaginary  parts  of  the  field.  Substituting  Eq.  (34)  into  Eq.  (33)  we 
obtain 

hi2)] 

:0  7?r  ''I2)  >2k  [«>*k2>  Id], 


OO 

x  Z 

k  = 


(35) 
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Curves  of  this  distribution  are  given  by  Beckmann  in  References  9  and  10.  The 
probability  density  function  of  the  power  can  be  obtained  if  desired  from  the  prob¬ 
ability  density  function  of  the  field  amplitude  by  the  relations  in  Reference  6 


P  |2  ( I  F  I  2) 

I  FT 


1 

IfF 


(36) 


using  the  notation  p^(x)  to  mean  the  probability  density  of  the  random  variable  y  as 
a  function  of  x. 


The  combinations  of  ^  and  —  appearing  in  Eqs.  (28)  and  (29)  for  the 
variance  of  the  real  and  imaginary  parts  of  the  field  can  be  replaced  by  their 
small  angle  forms  for  N.  >  4  thus  yielding 


aX*$T  [L  an  '  Z  an  cos  2(dnu+  %>] 

ay2  ST  [l  V  Z  3^c°s2(d  u+o  )  . 

L  n  n 


(37) 


(38) 


The  second  term  within  the  square  brackets  is  (for  real  fields)  the  function  H  <u' 

o 

encountered  earlier  (see  Eq.  (10)]  in  the  course  of  obtaining  an  expression  for  the 
variance  of  the  power.  We  showed  above  (see  p.  19)  that  H  (u)  can  be  neeleeted 

r-.  2  ° 

compared  to  unless  u  is  close  to  the  points,  midway  between  grating  lobes, 

for  which  |  H  (u)l  lakes  on  its  maximum  value.  Then 


9.  Beckmann,  P.  ,  and  Spizzichino,  A.  (1963)  The  Scattering  of  Electromagnetic 
Waves  Prom  Hough  Surfaces,  Pergamon  Press,  Macmillan,  New  York, 
Appendix  F. 

10.  Beckmann,  P.  (1963)  Statistical  distribution  of  the  amplitude  and  phase  of  a 
multiply  scattered  field,  ,1.  Res.  Natl.  Bur.  Std.  ,  D.  66(No.  31:231-240. 


3  1 


(39) 


ifi 


L2  2 

a  =  a 
n 


and  K=  1.  Since  Iq(0)=  1  and  1^(0)=  0  for  k^  0,  Eq.  (35),  which  gives  the  probability 
density  of  the  normalized  field  amplitude,  reduces  to 


p(  I  f  I  >  =  2  |  f  |  exp  -<B2  +  |f|2) 


.  (^4- 


p(  |  F| )  =  -U/-  exp 
a 


a  distribution  known  as  the  Rice-Nakagami  distribution  and  treated  in  slightly 
different  forms  in  detail  by  Rice  and  by  Norton  et  al  in  Reference  12. 


11.  Rice,  S.  O.  (1954)  Mathematical  analysis  of  random  noise,  reprinted  in 

Selected  Papers  on  Noise  and  Stochastic  Processes,  N.  Wax.  (Eds.  ), 
Dover,  New  York,  pp.  239-241. 

12.  Norton,  K.A.,  et  al  (1955)  The  probability  distribution  of  the  amplitude  of  a 

constant  vector  plus  a  Rayleigh-distributed  vector,  Proc.  IRE  43:  1354-1361. 

Norton  et  al  discuss  the  complement  of  the  cumulative  distribution  function  of  the 
probability  density  function 

P(D  =  ^  r  exp  [  -  (1  +  r2)/k2J  lQ  (  ~z) 
which  is  obtained  from  Eq.  (40)  via  the  correspondence 

I  i •  I  I  v  I  i,  i  2  o-2 


r  =  I  F 


|X|  +1,  k  =  2<j‘ 


Figure  4  of  Norton  et  al  contains  plots  of  the  complement  of  the  cumulative  distribu¬ 
tion  function  of  —  for  various  values  of  the  parameter  K=  20  log^k).  In  our 

notation  this  is  equivalent  to  plotting  the  complement  of  the  cumulative  distribution 

I  K  I  I  ^  I  2 

function  of  J  — —  for  various  values  of  K  =  10  log,n(2cr  ).  Rice  discusses  the 

\2  (J  10 

probabilitv  density  function 

'/  .2.2, 


p(v)  =  v  exp  I  -  - — ^ - ■  J  IQ  (av) 

which  can  be  obtained  from  Eq.  (40)  by  letting 


He  gives  plots  of  both  p(v)  and  its  cumulative  distribution  function. 


For  our  purposes  it  is  somewhat  more  convenient  to  follow  Rice's  normalization 


„M.«,  Jfi 

rr  *  n 


(41) 


For  small  values  of  the  parameter  a  the  cumulative  distribution  function 
v 


C(v)  v  exp  (-  — — j  IQ  (av) 


dv 


(42) 


can  be  obtained  by  numerical  integration.  For  large  values  of  a  Rice  obtains  an 
asymptotic  expansion  for  C(v)  in  inverse  powers  of  a  by  first  replacing  I^iav)  by 
its  large  argument  asymptotic  form,  given  in  Reference  8, 


a  v 


I  (av)  ~  "  — 

°  yrFav 


(l  +  av) 


so  that 


p(v)  ~ 


2  n 


nd  then  expanding^-  j  ^1+_av)  *n  a  ser*es  P°wers  of  and  integrating 


term  by  term.  The  resulting  expansion  through  terms  a  '  is 


C(v) 


2  ay[2~H 


1  - 


v  -  a 
4  a 


1  +  (v  -  a) 

iTo2 


exo\-{w-a)2  ' 
exp  2 


(43) 


In  Table  1  we  compare  values  of  the  Rician  cumulative  distribution  function  obtained 
by  performing  the  integration  in  Fq.  (42)  numerically  (using  the  IMSL  routine 
DCADRE)  with  values  obtained  using  the  asymptotic  form  of  Fq.  (43).  It  is  seen 
that  the  asymptotic  form  gives  an  excellent  fit  to  the  exact  values  for  a  >  3. 

If  we  further  specialize  our  results  and  consider  the  probability  distribution  of 
the  amplitude  at  a  direction  for  which  the  error-free  pattern  has  a  null,  then>v 
given  by  Fq.  (26)  equals  zero  and  p(  |  F  |  )  in  Fq.  (40)  reduces  to  the  Rayleigh 
distribution 


p( !  F  I  ) 


I  l'l 

— exp 

<7 
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According  to  Eq.  (36),  the  probability  density  function  of  the  power  is  then  the 
exponential  distribution 


The  mean  power  is  given  by 


=  -Vf 

2a~  Jo 


I-' I  2  e  2”  d([F|2)  =  202  =  A 


Comparing  Eq.  (46)  with  the  expression  in  Eq.  (14;  for  the  variance  of  the  power  at 
a  null  of  the  error-free  pattern  obtained  earlier  and  neglecting  Ho<u)in  Eq.  (14) 
to  be  consistent  with  the  assumption  made  in  obtaining  Eqs.  (39)  and  (40),  we  see 
that  Eq.  ( 14)  contains  the  extra  term  |  £  |wj4.  This  apparent  discrepancy  is 
explained  bv  the  fact  that  our  derivation  of  the  probability  distribution  of  the  ampli¬ 
tude  and  power  is  based  on  the  application  of  the  Central  Limit  Theorem  and  so  is 
an  approximation  to  the  actual  probability  distribution.  This  approximate  distribu¬ 
tion  approaches  the  actual  distribution  as  the  number  of  elements  in  the  array,  N, 
becomes  large,  provided  no  one  amplitude  of  the  element  weights  can  be  singled  out 
as  dominant.  Our  derivation  of  the  expression  for  the  variance  of  the  power 


(Eq.  (14)]  involved  no  approximation  and  so  cannot  be  expected  to  be  identical  to  the 

expression  for  the  variance  obtained  by  applying  the  Central  Limit  Theorem  except 

under  the  same  conditions  that  the  Central  Limit  Theorem  applies.  But  we  have 

seen  earlier  (seep.  20)  that  Y'  |w  I4  is  of  the  order  of  1/N  compared  with 
9  2 

(  V"!  w  |  )  for  any  amplitude  taper  described  by  a  polynomial  (that  is,  any  usual 
amplitude  taper)  and  indeed  for  such  tapers  no  one  amplitude  is  dominant.  Hence 
§  X]  I  wnl  4  can  be  neglected  compared  to  (  23 1  2)  under  exactly  the  conditions 
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for  which  the  Central  Limit  Theorem  applies,  thus  yielding  agreement  between 
Eqs.  (46)  and  (14).  Note  that  the  approximate  equality  of  the  variance  of  the  power 
and  the  square  of  the  mean  power  demonstrated  above  in  connection  with  Eq.  (16) 
approaches  actual  equality  as  the  number  ol  elements  in  the  array  increases  and  the 
probability  distribution  of  the  power  approaches  the  exponential  distribution  for 
which,  as  we  have  seen,  the  variance  is  exactly  equal  to  the  square  of  the  mean. 

When  nulls  are  placed  in  array  patterns  adaptively,  the  presence  of  multiple 
sources  of  interference  can  make  it  important  to  have  nulls  at  each  of  several 
different  directions.  Accordingly  it  is  of  interest  to  consider  the  effect  of  quantized 
phase  errors  on  the  pattern  at  several  directions  simultaneously,  for  each  of  which 
the  error-free  pattern  has  a  theoretically  perfect  null.  It  is  then  natural  to  examine 
the  distribution  of  the  least  deep  null  among  the  several  different  directions,  since; 
the  least  deep  null  may  serve  as  a  limiting  factor  in  determining  the  performance 
of  the  array.  Now,  tne  probability  that  the  least  deep  null  is  less  than  or  equal  to 
a  given  level  is  equal  to  the  probability  that  the  pattern  at  all  the  specified  direc¬ 
tions  is  less  than  or  equal  to  that  given  level.  If  we  make  the  simplifying  assump¬ 
tion  that,  for  a  given  distribution  of  phase  errors  across  the  array,  the  null  depths 
at  any  two  locations  for  which  the  error-free  pattern  has  a  null  are  independent  of 
each  other,  then  we  can  obtain  the  desired  probabilitv  simply  by  multiplying  together 
the  probabilities  that  each  null  is  less  than  or  equal  to  the  given  level.  Hut  from 
the  discussion  above  |see  Eq.  (44)],  the  probability  distribution  of  the  power  at  any 


locution  for  which  the  error  free  pattern  has  a  null  is  given  by 

I  |  ')  |  i 

a  i 3  -  Ail  -  f!~E 


ail-!3.  -  '  ( 

20  Jo 


d(!  !•':“)  =  1  -  e 


Hence  the  cumulative  probability  distribution  of  the  least  deep  null  among  A I  loca¬ 
tions  for  which  the  error-free  pattern  has  perfect  nulls  is  given  hv 

If  !2  x  .. 


Si(|1-'n '  = 


The  corresponding  probability  density  function  is 
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For  the  first  few  values  of  M  we  have  the  results  in  Table  2, 


Table  2.  Mean  and  Variance  of  the  Probability  Density  Function 
of  the  Least  Deep  Null  Among  M  Locations 


M 

S 

Var(?) 

1 

1 

1 

|f|2  =  2a2  T 

2 

3/2 

5/4 

2 

Var  ( |  F  |  2)  =  (2a2)  Var(0 

3 

11/6 

49/36 

2o2  =  7>2 

3  n 

4 

25/12 

205/144 

Note  that  for  M=  1,  the  values  for  the  mean  and  variance  of  the  least  deep  null  agree 
with  those  found  above  [compare  with  Eqs.  (45)  and  (46)]  for  a  single  null. 

The  final  special  case  we  shall  consider  here  is  that  of  the  probability  distribu¬ 
tion  of  the  amplitude  and  power  at  directions  midway  between  grating  lobes.  For 
a  linear  phase  variation,  <j>n  =  -dnug,  of  the  element  weights  across  the  arrav,  these 
directions  are  given  by 

u  =  u  ±  m  n  . 
s 

Since  the  focus  of  this  report  is  the  influence  of  phase  errors  on  sidelobe  levels  we 
will  restrict  our  attention  to  the  case  of  m  odd.  (If  m  is  even  there  is  a  grating 
lobe  at  the  location  u.  The  assumption  that  m  is  odd  is  not  immediately  needed 
here  but  is  required  in  demonstrating  agreement  between  the  exact  expression  for 
the  variance  of  the  power  derived  earlier  and  the  expression  derived  below  from 
the  probability  distribution  of  the  amplitude  and  power  for  this  special  case.  )  For 
these  directions  it  is  easily  verified  that 

y  a2  cos  I  2  (0  +d  u)  =  ±  y  a2  (49) 

nn  LnnJ  nn 


so  that  from  the  approximations  (37)  and  (38),  either 
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n 


(50) 


(51) 


(j  r  =  o 
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according  as  the  plus  or  minus  sign  respectively  is  taken  in  Eq.  (49).  Rather  than 
attempting  to  specialize  the  general  expressions  we  have  derived  above  for  the 
probability  distribution  of  the  amplitude  and  power  to  cover  these  cases,  it  is  much 
simpler  to  derive  the  desired  distributions  starting  with  the  variances  of  the  real 
and  imaginary  parts  of  the  field,  Eqs.  (50)  or  (51).  (The  means  of  the  real  and 
imaginary  parts  of  the  field  are,  from  Eqs.  (26)  and  (27),  given  by  — Fq(u) 
and  0  respectively.  ) 

For  the  first  case,  Eq.  (50),  the  amplitude  of  the  field  is  given  by 
|  F  |  =  (X2  +  Y  “) 1  ^2  so  that  Y  =  ±  ( |  F  |  2  -  X2)  ^  ^  2.  Using  the  formula  for  the  probability 
distribution  of  a  function  of  a  random  variable^  and  the  fact  that  Y  is  Gaussian 
distributed  we  have  immediately 


I  2  -  >T2 


P  ( I  F  |  )  = 


_1_  /  2\» 

%  {  *  ) 


( I  F  I  -  X  > 


The  probability  density  of  the  power  is.  from  Eq.  (36),  given  by 


I  !•’  I  2  -  X  2 


Vi 2 


(If  |2)  = 


f|2  2  x2 


(2„)w  CTV  (|Fr-X‘)‘ 


and  a  short  calculation  then  gives  the  cumulative  distribution  function  of  the  power, 
r / ,  . o 


If!2  >  X2 


where  erf(z)  is  the  error  function. 


erf(z)  =  -L  f 
■JT  J o 


e  *  d  t 


The  mean  power  is  given  by 


( I  K I  2  -  X2) 


FFF  \n  ~  ^ 
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Making  the  substitution  v  =  |  F  |  2  -  X2,  dv  =  d  ( |  F  |  2).  and  using  the  formulas  given 
in  Reference  13 


in  agreement  with  Eq.  (5). 

The  variance  of  the  power  is  given  by 


If)  2  -x2 
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2a, 
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2  (|l-i 


2-x2)T7T 


d  ( |  F  |  2  >  -  (X2+ay2)2 


Making  tne  same  cnange  of  variables  used  in  obtaining  the  mean  and  again  using 
Eq.  (53)  we  find  tnat  tne  variance  of  tne  power  is 

Var(|f'|2)  =  2ay4 


Comparing  t lie  approximation  (55)  with  tne  exact  (to  within  tue  small  phase  error 

approximation)  expression  for  th’e  variance  of  tne  power  given  by  Eq.  (13)  and  noting 

that  for  the  special  case  we  are  considering  H  (u)  =  51  a2  cos  ( 2  d  (u  -  u  )]  =  51  a  , 

1  “o  n  ns  n 

we  see  that  tne  exact  expression  contains  several  terms  in  addition  to  those  in  the 

approximation  (55).  Similarly  to  our  discussion  on  p.  35  of  tne  apparent  discrepancy 

at  a  pattern  null  between  tne  exact  expression  for  the  variance,  and  the  expression 

obtained  from  tne  Central  Limit  Theorem  and  the  probability  distribution  of  the 

amplitude,  it  can  be  snown  (see  Appendix  A)  that  the  extra  terms  can  be  neglected 

for  large  arrays  and  so  the  two  results  agree. 

For  the  second  case,  we  consider  the  approximation  (51),  and  begin  by  noting 

tnat  the  error-free  pattern  has  a  null  at  the  directions  we  are  considering.  Since 
2  2 

51  a  cos  (2d  (u-u  )]=  -  T"  a,  cos  (2d  <u  -  u  )]  =  - 1  and  hence 
^  n  1  n  s1  ^  n  'n  ?  ‘ 


40 


cos  [2d  (u  -  u  )]  =  cos 
1  n  s 


1 2  -  (n-  l)j  mtlj 


=  cos  [(N  -  1)  m  v  ] 


so  that 


(N  -  l)mi  =  ±  n,  ±  3u,  ±  5tt, 


(N  -  l)m  =  ±  1,  ±3.  *  5, 


But  then 


cos  Id  (u  -  u  )]  =  cos 
1  n  s 


|  -  (n  -  dJ  m  77  J 

=  (-  D<n  '  1)m  Cos(^  m*) 


since  (N  -  D  m  is  odd  and  so 


F  (u)  =  /J  a  cos  |d  (u  -  u  )]  =  0 
o  n  n  s' 


Since  the  error-free  pattern  has  a  null,  the  mean  of  the  real  part  of  t lie  field, 
”K,  equals  zero.  Thus  |  F  I  =  |x|  and  the  probability  density  of  the  amplitude  and 


power  are 
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respectively.  These  density  functions  are  identical  with  those  obtained  for  the  first 


case  bv  the  substitutions  X 
of  the  power  are 
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respectively.  Comparing  the  approximation  (56)  with  the  more  exact  expression 
given  by  Eq.  (13)  for  the  variance,  noting  that  all  the  terms  in  the  approximation 
(13)  containing  Fo(u)  are  zero  because  of  the  null  in  the  error  free  pattern  and  that 

I  Hq(u)  I  2  a2  j  ,  and  recalling  that  ^  a4  is  of  the  order  of  1/N  compared  to 

^^a2  j  for  large  N,  we  see  that  the  two  expressions  for  the  variance  are  in 


agreement . 


3.  NUMERICAL  RESULTS  AND  DISCUSSION 

Example  calculations  in  this  section  give  results  that  are  compared  with  re¬ 
sults  obtained  in  the  analysis  of  the  preceding  section.  We  consider  here  the  dis¬ 
tribution  of  the  power  at  (1)  a  single  null,  (2)  multiple  nulls,  and  (3)  a  non-null 
sidelobe  location. 

3. 1  Single  Null 

We  begin  with  the  distribution  of  the  power  at  a  direction  for  which  the  error- 
free  pattern  has  a  perfect  null.  For  concreteness  we  consider  a  broadside  array  of 
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1 

'  ■*- 


•  1 


1 
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79  elements  with  half  wavelength  spacing  and  a  40  dB  Chebyshev  taper.  The  phases 
of  the  array  elements  are  assumed  to  be  set  with  8  bit  digital  phase  shifters.  The 
error-free  pattern  of  the  array  has  a  null  at  20.  3989. .  .  “  between  the  13th  and  14th 
sidelobes.  Calculating  the  mean  power  at  this  location  using  Eq.  (4)  we  find 

T7P  =  0.  8072  X  10'8  =  -60.  9  dB 

and  calculating  the  variance  of  the  power  from  Eq.  (12), 

Var  ( |  F  J  2)  =  0.  6351  X  10~12  ; 


hence,  the  standard  deviation  of  the  power  equals  0.7969  X  10  .  Note  that  the 

standard  deviation  of  the  power  is  approximately  equal  to  the  mean  power,  in  agree¬ 
ment  with  the  discussion  following  Eq.  (16).  The  contribution  to  the  variance  of 
first  three  terms  on  the  right  hand  side  of  Eq.  (12)  is  effectively  ze ro  because  of 

the  null  of  the  error-free  pattern;  the  contribution  of  the  fourth  term  (in  |h  (u)|2) 

-17  -12  ° 

equals  0.  84  X  10  ,  the  fifth  term  equals  0.  65  X  10  ,  and  the  last  term  equals 

0.  17  X  10  18.  Thus,  as  was  shown  on  p.  20,  |h  (u)|2  is  negligible  compared  to 

/  V  l  i  2  \2 

l  /Jw  I  (except  at  directions  halfway  between  grating  lobes),  and,  as  dis- 

'  n  /  |  1 4 

cussed  on  p.  20  and  proved  in  Appendix  A,  the  term  in  7,1  w  |  is  of  the  order  of 

2  n 

l/N  compared  with  the  term  in  (y  )  • 

We  compute  the  probability  distribution  of  the  power  at  the  same  location  from 

Eqs.  (28)  and  (29)  with  the  result,  CT^.  =  0.  405  X  10  8  and  =  0.  402  X  10  8  with 

the  term  in  T"]  a2  equal  to  0.403  X  10"®  and  the  term  in^a2  cos  {2x1/  )  equal  to 
_ ^  n  n  2  n 

0.  145  X  10  and  hence  negligible  compared  to  the  term  ^  a^  as  was  shown  in  the 

discussion  of  Eq.  (38).  The  parameter  B  of  the  normalized  probability  distribution, 

Eq.  (34a),  is  effectively  zero  because  of  the  pattern  null,  and  the  parameter  K  given 

bv  Eq.  (34c)  is  equal  to  0.  996.  Since  K  =  1  and  X  =  0,  we  expect  the  probability 

distribution  of  the  power  to  be  well  approximated  bv  the  exponential  distribution 

A 2  2  -6 

given  in  Eq.  (44)  with  a  mean  of  — j-  2^an  ~  X  10  and  a  variance  equal  to 

the  square  of  the  mean,  0.  65  16X  10"  Note  that  the  variance  of  the  power  cal¬ 
culated  assuming  the  exponential  distribution  of  the  power  is  quite  close  to  the  value 
of  0.  6351  X  10  6  obtained  from  the  exact  expression  for  the  variance. 

To  verify  the  above  theoretical  results  we  wrote  a  computer  program  to  ran¬ 
domly  perturb  the  phases  of  an  equispaced  array  and  calculate  the  resulting  power 
at  a  specified  angle,  given  the  number  of  elements,  the  spacing,  the  amplitude 
taper,  and  the  number  of  bits  of  the  phase  shifters.  Each  phase  is  independently 

perturbed  from  its  error-free  value  bv  a  random  number  taken  from  a  uniform 

N  \ 

distribution  in  the  interval  (-»/ 2  ^>’t.  nl 2  A  sample  size  nf  1000  patterns  was 


I 


taken  for  a  79  element,  half  wavelength  spacing,  40  dB  Chebyshev  tapered  array. 

The  mean  power  at  the  location  6  =  20.  3989.  . .  °  was  found  to  be  0.  8043  X  10  ^ 

-  1 2 

and  the  variance  to  be  0.  6167  X  10  ,  in  close  agreement  with  the  theoretical  re¬ 

sults.  In  figure  1  we  show  a  computer  generated  plot  (obtained  using  the  IMSL 
routine  t.'SPC)  of  the  sample  cumulative  probability  distribution  and  the  theoretical 
distribution  of  Eq.  (47).  (The  letter  M  indicates  that  both  curves  share  the  same 
print  location.  )  The  close  fit  of  the  two  curves  is  apparent.  To  further  check  the 
agreement  of  the  sample  and  theoretical  cumulative  distributions,  a  Kolmogorov- 
Smirnov  test  (described  in  Reference  14)  was  performed  using  the  IMSL  routine 
NKS1.  The  maximum  absolute  difference  between  the  sample  cumulative  distribu¬ 
tion  and  the  theoretical  exponential  distribution,  (n=  size  of  sample,  here  equal 
to  1000)  was  found  to  be  0.  0228,  the  statistic  z  =  V*  D  =0.7  2,  and  the  probability 
of  obtaining  a  value  of  z  equal  to  or  greater  than  0.  72  was  found  to  be  0.  68  so  that 
there  is  no  reason  to  reject  the  hypothesis  that  the  exponential  distribution  is  the 
underlying  population  distribution  of  the  power. 

A  similar  set  of  calculations  was  performed  for  a  nu,!  that  was  imposed  on  the 
79  element,  40  dll  Chebyshev  pattern  by  small  phase  perturbations,  using  an  itera¬ 
tive  technique  described  in  Reference  15.  The  results  do  not  differ  in  any  significant 
way  from  those  obtained  for  the  "natural"  null  and  need  not  be  described  here. 

for  reference  purposes  it  is  useful  to  have  a  plot  of  the  normalized  exponential 
cumulative  probability  distribution  of  the  power  at  a  direction  for  which  the  error- 
free  pattern  has  a  null.  In  figure  2  we  have  plotted  the  cumulative  probability 
distributions  [  Eq.  (48)] 


CM<p2*  =  ^1  -  exp  (-  p2)  j 

2  2 

for  M  =  1  to  3  as  a  function  of  10  log  n  (p  )  where  p  is  the  power  normalized  by 

o  1 U 

2  a  with 


The  curve  M  =  1  is  the  exponential  cumulative  probability  distribution.  (The  curves 
for  other  values  of  M  are  relevant  to  the  discussion  below  of  multiple  nulls.  )  As 
an  example  of  the  use  of  this  plot,  in  the  case  considered  above  N,  .,=  8  so  that 


(Because  of  the  length  of  Reference  14,  References  14  and  15  will  not  be  listed  here. 
Sec  References,  page  57.  ) 


A=  ^ / 256  51  The  sum  norma^zec^  by(y"]an  j  is  calculated  to  be 

0.  01608  so  that  2(7^  =  0.  807  2  X  10  ^  =  -60,  9  dB.  The  probability  of  a  null  depth 
2 

20  dB  below  2a  (-80.  9  dB  from  the  mainlobe  peak)  is  then  seen  to  be  0.  01,  and 

2 

the  probability  of  a  null  depth  10  dB  below  2a  (-70,  9  dB  from  the  mainlobe)  is 

2 

0.  095.  The  probability  of  a  value  of  power  more  than  -4  dB  below  2a  (that  is, 
more  than  -64.  9  from  the  mainlobe)  is  0.  5. 
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Figure  1.  Sample  (1)  and  Theoretical  (2)  Cumulative  Distribution  Functions  of 
Power  at  the  Error-free  Pattern  Null  Location  9  -  20.  3989.  . .  °.  Half  wave¬ 
length  spacing  is  used  in  a  1 9-element,  40  dB  Chebyshev  tapered  array  with 

phase  errors  uniformly  distributed  in  |-A.  A],  A  =  =  0.012  (M  indicates 

that  both  curves  share  the  same  print  locution) 
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Figure  2.  The  Cumulative 
Probability  Distribution  Function 

Cw(p2)  =  |  1  -  exp  (-  p2)]W  for 
M  =  1,  2,  and  3 


3.2  Multiple  Nulls 

We  next  describe  computations  to  compare  with  the  analytic  results  obtained 
for  the  distribution  of  the  least  deep  null  among  several  directions  for  each  of 
which  the  error-free  pattern  has  a  perfe’t  null.  We  showed  in  the  previous  section 
that,  if  we  make  the  simplifying  assumption  that  the  null  depths  at  any  two  different 
locations  are  statistically  independent  of  each  other,  the  cumulative  distribution  of 
the  least  deep  null  among  M  locations  is  given  by  Fq.  (48).  In  the  first  computa¬ 
tion  performed  here,  we  took  a  sample  of  1000  random  perturbations  of  the  phases 
of  the  79  element,  half  wavelength,  40  dB  Chebyshev  array  with  N^t  =  8,  and 
examined  the  distribution  of  the  least  deep  null  between  the  locations  9  ~  20.  3989.  .  .  0 
and  21.  9620.  . .  °.  These  locations  are  the  nulls  between  the  13th  and  14th  sidelobes, 
and  between  the  14th  and  15th  sidelobes,  respectively.  The  values  predicted  for 
the  mean  and  variance  of  the  sample  on  the  basis  of  Fq.  (48)  (see  Table  2)  are 

|  F | 2  =  J  (2a2)  =  2  (0.  8072  X  10'6) =  0.  121  1  v  10* 5 


and 
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Var  ( j  F  |  2)  =  |  (2a2)  =  0.  8145  X  10‘12  . 

The  values  calculated  from  the  sample  of  1000  were 

|f|2100Q  =  0.  1123  X  10'5.  Var  <1^1  2>1000  =  0.7297  X  10'12  , 

reasonably  close  to  the  predicted  values  although  differing  somewhat  more  in  rela¬ 
tive  terms  than  the  corresponding  results  for  a  single  null  at  8  =  20.  3989.  . .  0 
described  above.  The  computer  generated  plot  of  the  sample  and  theoretical 
cumulative  distribution  functions  is  shown  in  Figure  3.  Comparing  Figure  3  with 
Figure  1  shows  that  the  fit  of  the  theoretical  with  the  sample  cumulative  distribu¬ 
tion  function  (CDF)  is  somewhat  poorer  for  the  two-null  case  especially  in  the 
middle  portion  of  the  plot.  Note  that  the  sample  CDF  is  fairly  consistently  higher 
than  the  theoretical  CDF  indicating  that  a  greater  proportion  of  the  sample  values 
of  the  least  deep  null  are  to  be  found  in  a  lower  part  of  the  range  of  power  than  is 
predicted  by  the  theoretical  distribution.  A  Kolmogorov-Smirnov  test  was  per¬ 
formed  on  the  sample.  The  maximum  absolute  difference  between  the  sample  and 
theoretical  CDF's,  D^,  was  0.  07,  the  statistic  z  =  2.  21,  and  the  probability  of 
obtaining  a  value  of  z  equal  to  or  greater  than  2.  21  was  found  to  be  0.  0001  indicating 
that  the  hypothesis  that  the  underlying  population  distribution  of  the  sample  is  des¬ 
cribed  by  the  theoretical  distribution  given  by  Eq.  (48),  should  be  rejected. 

A  similar  computation  was  performed  for  the  case  of  two  nulls  imposed  on  the 
7  9  element,  40  dB  Chebyshev  pattern  by  small  phase  perturbations  at  the  locations 
6  =  20.  78°  and  21.  2°  (approximately  the  left  -3  dB  point  and  the  sidelobe  peak  of 
the  13th  sidelobe).  The  sample  mean  and  variance  were 

|F|21000  =  0.9044  X  10~6,  Var(|F|2)l0()0=  0.7196  X  10'U. 

A  plot  of  the  sample  and  theoretical  CDF's  is  shown  in  Figure  4.  Here  the  diver¬ 
gence  between  the  two  distributions  is  clearly  apparent,  the  sample  distribution 
being  significantly  higher  than  the  theoretical  distribution  over  almost  the  entire 
range.  This  of  course  is  reflected  in  the  considerably  lower  value  of  the  sample 
mean  compared  with  the  theoretical  mean.  The  Kolmogorov-Smirnov  test  gave 

D  =  0.  2080,  z  =  6.  58,  and  Prob  z  2  6.  58  =  0.  0000. 
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Figure  3.  Sample  (1)  and  Theoretical  (2)  Cumulative  Distribution  Functions  of 
the  Maximum  Value  of  Power  at  the  Two  Error-free  Pattern  Null  Locations 
9  =  20.  3989. .  .  °  and  21.  9620. .  .  °.  Half  wavelength  spacing  is  used  in  a 
79-element  40  dB  Chebyshev  tapered  array  with  phase  errors  uniformly 

distributed  in  [-A,  A],  A  =  =  0.  012  (M  indicates  that  both  curves  share 

the  same  print  location) 


The  fact  that  the  sample  CDF  is  consistently  higher  than  the  theoretical  CDF 
indicates,  as  mentioned  above,  that  a  significantly  greater  proportion  of  the  sample 
values  for  the  least  deep  null  are  below  those  predicted  by  the  theoretical  distribu¬ 
tion,  and  suggests  that  the  null  depths  at  closely  spaced  locations  of  the  pattern  are 
not  independent  of  each  other  as  was  assumed  in  deriving  the  theoretical  CDF"  given 
by  Eq.  (48).  If  the  null  depth  is  low  at  one  location  it  tends  to  be  low  at  the  other 
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Figure  4.  Sample  (1)  and  Theoretical  (2)  Cumulative  Distribution  Functions  of 
the  Maximum  Value  of  Power  at  the  Two  Error-free  Pattern  Imposed  Null 
Locations  9  —  20.  78°  and  21.  2°.  Half  wavelength  spacing  is  used  in  a 
79-element,  40  dB  C'hebyshev  tapered  array  with  phase  errors  uniformly 

distributed  in  [-A,  A],  A  =  =  0.  012  (M  indicates  that  both  curves  share 

the  same  print  location)  430 


location  as  well.  As  a  check  on  the  validity  of  this  hypothesis,  another  calculation 
was  performed  in  which  nulls  were  imposed  on  the  basic  7  9  element,  40  dB 
C'hebyshev  pattern  at  the  widely  spaced  locations  6  =  12°  and  72°,  and  a  sample 
taken  of  1000  values  of  the  least  deep  null  among  these  two  locations.  The  sample 
mean  and  variance  were 
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!  l'i  ~l000  =  0.  1  177  ■  10-  (-59.3  dB),  Var(|F|2)  =  0.7747  X  10-12  , 

and  the  plot  of  the  sample  and  theoretical  CDF's  is  shown  in  Figure  5.  The  much 
closer  fit  of  the  sample  and  theoretical  CDF's  than  in  Figure  4  is  apparent  and  is 
reflected  in  the  Kolmogorov-Smirnov  test  which  gave 

D  =  0.  025.  z  =  0.  776,  Prob  z  2  0.  7  76  =  0.  58. 
n 

indicating  consistency  of  the  sample  and  theoretical  CDF's. 
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Figure  5.  Sample  (1)  and  Theoretical  <2)  Cumulative  Distribution  F  unctions  of 
the  Maximum  Value  of  Power  at  the  Two  Krror-free  Pattern  Imposed  Null 
Locations  9  =  12°  and  72°.  Half  wavelength  spacing  is  used  in  a  79-element, 
40  dH  Chebyshev  tapered  array  with  phase  errors  uniformly  distributed  in 

(-A,  A|,  A  =  =  o.  012  (M  indicates  that  both  curves  share  the  same  print 

location)  b 
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The  high  statistical  correlation  between  the  null  depths  at  neighboring  loca¬ 
tions  has  interesting  implications  for  nulling  with  phased  arrays  employing  digital 
phase  shifters  when  closely  spaced  nulls  are  desired.  It  means  that  the  assumption 
of  independence  of  null  depths  gives  an  overly  conservative  estimate  of  the  prob¬ 
ability  that  the  null  depths  at  two  or  more  closely  spaced  locations  (for  each  of 
which  the  error-free  pattern  has  very  deep  nulls)  will  both  be  less  than  a  given 

level.  Suppose,  for  example,  that  the  probability  of  a  single  null  being  20  dB  below 
2 

2cr  is  0.  1.  Then  under  the  assumption  of  independence  we  would  conclude  that  the 

2 

probability  that  two  adjacent  nulls  were  both  less  than  20  dB  below  2a  was  0.  01 

(see  Figure  2  in  which  t lie  theoretical  CDF,  Eq.  (48),  is  plotted  for  n=  1,  2,  and  3). 

In  fact,  however,  things  are  not  all  that  bad  and  the  probability  that  both  nulls  are 

2 

less  than  20  dB  below  2 a  may  well  be  closer  to  0.  I  than  0.  01.  To  illustrate  the 
fact  that  the  one  null  theoretical  CDF  may  give  a  better  fit  to  the  two  null  sample 
data  than  the  two  null  theoretical  CDF,  in  Figure  6  we  plot  the  sample  CDF  of  Fig¬ 
ure  4  (corresponding  to  nulls  at  20.  78°  and  21.  2°)  along  with  the  single  null 
theoretical  CDF.  The  single  null  theoretical  CDF  is  clearly  closer  to  the  sample 
CDF  than  is  the  two-null  CDF  (z=  2.  65  under  the  assumption  of  the  one-null  CDF 
whereas  z=  6.  58  for  the  two-null  CDF).  We  also  performed  a  similar  computa¬ 
tion  for  the  case  of  three  nulls  at  9  =  20.  78°,  21.2°,  and  21.8°,  and  found  that  the 
one -null  theoretical  CDF  gave  a  better  fit  to  the  three -null  sample  distribution  than 
either  the  three-null  or  two-null  theoretical  CDF's  (z=  4.  59  under  the  assumption 
of  the  one-null  theoretical  CDF,  whereas  z=  9.  43  for  the  three-null  CDF  and  4.  89 
for  the  two-null  CDF).  It  would  be  interesting  to  study  in  detail  the  quantitative 
dependence  of  the  correlation  of  null  depths  on  the  angular  separation  of  the  null 
locations  and  on  other  factors  such  as  the  magnitude  of  the  quantization  error,  the 
number  and  spacing  of  the  elements,  and  the  amplitude  taper.  Such  an  investiga¬ 
tion,  however,  lies  beyond  the  scope  of  the  present  study  and  we  limit  ourselves 
lie  re  to  drawing  attention  to  the  fact  that  such  a  correlation  is  important  in  deter¬ 
mining  the  null  depths  in  the  error  pattern  at  closely  spaced  null  locations. 

3.3  Non-null  SidelolK*  location 

To  verify  the  results  obtained  in  the  previous  section  for  the  theoretical  dis¬ 
tribution  of  power  at  a  non-null  sidelobe  location  not  halfway  between  grating  lobes 
[see  the  discussion  of  Eq.  (38)],  a  sidelobe  location  of  8  =  20.  1°  was  chosen  in  the 
13th  sidelobe  of  the  7  9  element,  40  dB  Chebvshev  pattern.  In  Figure  7  we  show  the 
computer  generated  plot  of  the  sample  CDF  based  on  1000  samples  of  phase  errors 
with  =  8,  and  the  theoretical  Hician  CDF  calculated  using  the  asymptotic  form 
(43).  For  this  value  of  8,  the  parameter  a=  2o“X  =  8.  99  so  that  the  asymptotic  form 
is  an  excellent  approximation  to  the  exact  integral  of  the  Hician  probability 
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distribution  function  The  close  fit  of  the  sample  and  theoretical  distribution  func¬ 
tions  is  apparent  an,d  is  corroborated  by  the  results  of  the  Kolmogorov-Smirnc v 

test:  =0.02,  z=-0.  63.  and  Prob  z  >  0.  63  =  0.  82. 
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f  igure  6.  Sample  (1)  and  Theoretical  (2)  Cumulative  Distribution  Functions  of 
the  Maximum  Value  of  Power  at  the  Two  Error-free  Pattern  Imposed  Null 
Locations  6  =  20.  78°  and  21.  2°;  Theoretical  CDF  is  for  One  Error-free  Pattern 
Null  Location.  Half  wavelength  spacing  is  used  in  a  79-element,  -40  dli 
Chebyshcv  tapered  array  with  phase  errors  uniformly  distributed  in 

[-A,  A],  A  =  , , "  =  0.  012  ( A I  indicates  that  both  curves  share  the  same  print 

locations)  3  ’ 
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Figure  7.  Sample  (1)  and  Theoretical  (2)  Cumulative  Distribution  Functions  of 
Power  at  the  Lrror-free  Pattern  Non-null  Location  9  =  20.  1°.  Half  wavelength 
spacing  is  used  in  a  79-element,  40  dB  Chebyshev  tapered  array  with  phase 

errors  uniformly  distributed  in  (-A,  A],  A  =  =  0.012  (M  indicates  that  both 

curves  share  the  same  print  location) 


For  the  case  in  which  a  non-null  sidelobe  location  is  halfway  between  grating 
lobes,  the  theoretical  mean,  variance,  and  C'DF,  are  given  respectively  by  the 
approximations  (54)  and  (55),  and  bv  Fq.  (52).  To  compare  with  these  theoretical 
results,  a  sample  distribution  was  generated  starting  with  a  79  element,  one  wave¬ 
length  spacing,  40  dB  Chebvshev  taper  broadside  array,  and  introducing  random 

N  '  N... 

phases  uniformly  distributed  in  the  interval  |-tt/2  1  ,  77/2  ]  with  N.  .  =  8.  The 
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sample  mean  and  variance  of  the  power  at  the  locations  9  =  30°,  halfway  between 
grating  lobes,  based  on  a  sample  size  of  1000,  wereO.  1019  X  10  ^  and  0.  1345  X  10 
respectively  as  compared  with  the  theoretical  values  of  0.  1018  X  10  ^  (-39.  96  dB) 
and  0.  1303  X  10  **.  In  Figure  8  we  show  computer  generated  plots  of  the  sample 
and  theoretical  CDF's.  The  two  plots  appear  to  agree  well  but  this  apparent  agree¬ 
ment  is  not  corroborated  by  the  Kolmogcrov-Smirnov  test  which  gives  Dfl  =  0.  075, 
z=  2.  37,  and  Prob  z=5  2.  37=  0.  0000.  The  explanation  of  this  discrepancy  between 
the  result  of  the  Kolmogorov-Smirnov  test  and  the  close  agreement  of  the  theoretical 
and  sample  mean,  variance,  and  distribution  function  plots,  lies  in  the  fact  that  the 
theoretical  results  are  based  on  the  approximations  (50)  which  preclude  any  varia¬ 
tion  of  the  real  part  of  the  field.  As  a  result,  the  theoretical  power  distribution 
has  no  values  of  power  less  than 


In  fact,  however,  the  real  part  of  the  field  has  a  very  small  but  non-zero  variance 
(if  the  variance  of  the  real  and  imaginary  parts  of  the  field  are  calculated  using 
llqs.  (28)  and  (29)  we  find  that  =  0.  8104  X  10*11  and  o^.  =  0.  8072  X  10-6)  so  that 
a  number  of  the  sample  values  of  power  are  slightly  smaller  than  the  theoretical 
lower  limit.  This  results  in  a  significant  oifference  between  the  sample  and 
theoretical  distribution  functions  at  the  very  beginning  of  the  range  of  power  which, 
although  not  apparent  in  the  plots,  makes  itself  felt  as  a  sufficiently  large  value 
ol  It  to  result  in  a  Kolmogorov-Smirnov  rejection  of  the  hypothesis  that  the  sample 
ot  distribution  of  power  is  well  described  bv  the  theoretical  distribution.  Apart 
from  tins  difference  in  behavior  at  the  lower  limit  of  the  range  in  power,  the 
theoretical  and  sample  distributions  agree  well  as  is  evidenced  by  the  plots  and  by 
the  respective  values  for  the  mean  and  variance. 

t  <mrn  siov* 


In  this  retort  hc  have  analyzed  the  statistical  distribution  of  sidelobe  power  of 
a  linear  arrav  of  isotropic,  equispaced  elements  whose  excitation  coefficients  are 
subn  et  to  random  phase  errors.  Since  our  study  was  motivated  bv  the  desire  to 
describe  the  errors  introduced  in  the  array  pattern  by  the  use  of  digital  phase 
shifters  to  set  the  excitation  phases,  the  random  phase  errors  were  assumed 

Nfait 

throughout  to  be  uniformly  distributed  in  the  interval  [-n/2  ,  n/2  1  ]  where 

N,.,  is  the  number  of  bits  in  the  phase  shifters. 
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Figure  tf.  Sample  (1)  and  Theoretical  (2)  Cumulative  Distribution  Functions  of 
the  Power  at  the  Location  6  =  30°,  Midway  Between  Grating  Lobes  of  the 
Frror-free  Pattern.  One  wavelength  spacing  is  used  in  a  79-element,  40  dB 
Chebyshev  tapered  array  with  phase  errors  uniformly  distributed  in 

|-A,  A],  A  =  =  0.  012  (M  indicates  that  both  curves  share  the  same  print 

location) 
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We  began  by  deriving  general  expressions  for  the  mean  and  variance  of  the 
power  at  an  arbitrary  pattern  location,  and  obtained  small  phase  error  approxi¬ 
mations  valid  for  >4.  We  showed  that  for  large  arrays,  the  variance  of  the 
power  at  a  sidelobe  location  for  which  the  error-free  pattern  has  a  null  is  approxi¬ 
mately  equal  to  the  square  of  the  mean  of  the  power. 

We  then  derived  a  general  expression  for  the  probability  distribution  of  the 
amplitude  and  power  of  the  error  field  under  the  assumptions  that  the  error-free 
field  was  real  and  that  the  number  of  elements  in  the  array  was  sufficiently  large 
for  the  Central  Limit  Theorem  to  be  applied.  This  general  expression  was  then 
specialized  to  the  case  of  small  phase  errors  and  ordinary  sidelobe  locations  (loca¬ 
tions  not  halfway  between  grating  lobes).  For  this  case  the  distribution  of  sidelobe 
amplitude  is  described  by  the  Rice-Nakagami  distribution.  At  a  pattern  location 
for  which  the  error-free  pattern  has  a  null,  the  Rice-Nakagami  distribution  reduces 
to  the  Rayleigh  distribution,  and  the  probability  distribution  of  the  power  is  then 
described  by  the  exponential  distribution.  The  expressions  for  the  mean  and 
variance  of  the  power  obtained  using  the  exponential  distribution  were  compared 
with  the  general  expression  for  the  mean  and  variance  of  the  power  obtained  earlier, 
and  the  two  sets  of  expressions  shown  to  be  consistent  under  the  assumption  of 
large  arrays  used  in  deriving  the  probability  distribution  of  the  power.  Assuming 
statistical  independence  of  the  depth  of  null  at  two  or  more  sidelobe  locations  for 
which  the  error-free  pattern  has  a  null,  an  expression  for  the  statistical  distribu¬ 
tions  of  the  least  deep  null  among  several  locations  was  obtained  by  multiplying 
together  the  distributions  for  the  individual  nulls.  Finally,  the  distribution  of  power 
at  locations  halfway  between  grating  lobes  was  obtained  starting  with  expressions 
for  the  variance  of  the  real  and  imaginary  parts  of  the  field,  and  showing  that  at 
such  points  the  variance  of  either  the  real  or  imaginary  part  of  the  field  is  approxi¬ 
mately  equal  to  zero. 

Following  the  theoretical  analysis,  Monte  Carlo-type  computer  simulations 
were  performed  to  compare  with  the  theoretical  results.  Extensive  use  was  made 
of  the  Kolmogorov-Smirnov  test  to  test  the  hypothesis  that  the  sample  distribution 
of  power  generated  by  the  simulation  could  be  regarded  with  high  probability  as 
being  drawn  from  the  corresponding  underlying  distribution  obtained  theoretically. 

In  general,  close  agreement  between  the  simulation  and  theoretical  results  was 
obtained.  However,  the  assumption  of  statistical  independence  of  null  depths  at 
several  locations  for  which  the  error-free  pattern  has  a  null  was  shown  to  be 
invalid  for  closely  spaced  nulls.  The  correlation  between  the  null  depths  at  closely 
spaced  locations  increases  the  probability  of  the  shallowest  null  having  a  given 
minimum  depth  from  the  independent -multiple  null  probability  to  a  value  close  to 
the  probability  of  a  single  null  of  that  depth. 
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Appendix  A 

Proof  of  Theorems  Used  in  the  Analysis 


In  this  appendix  we  prove  two  results  used  in  the  analysis  section  of  the  report. 
The  first  result,  made  use  of  on  pp.  20  and  35,  is  that  for  any  taper  of  the 
amplitude  Iw^l  of  the  element  excitation  coefficients  describable  by  a  polynomial 
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n=  1 


~pt 


n-  1 


(Al) 


To  prove  (Al)  we  start  with  the  assumption  that  the  amplitudes  of  the  excitation 
coefficients  can  be  expressed  as  a  polynomial  of  degree  M, 


w  i  = 
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ak  (1 


2  n"  1  > 


n  =  1,2,, 
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(A2) 


with  M  and  the  a^  independent  of  N.  Eq.  (A 2)  in  effect  defines  a  family  of  element 

amplitudes,  each  member  of  the  family  corresponding  to  a  particular  value  of  N. 

The  form  of  the  argument  of  the  polynomial  is  chosen  to  equal  +1  and  -1  for  n=  1  and 

n- 1  k 

N  respectively.  Expanding  Eq.  (A2)  and  collecting  terms  in  ( tt — t)  we  can  write 
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as  claimed. 

The  second  result,  used  in  the  discussion  of  the  approximation  (55),  is  that  at 
a  pattern  direction  midway  between  grating  lobes. 


u  =  u  ±  m IT,  m  odd 
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for  which 
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H  (u)  =  57  a^  cos  [2d  (u  -  u  )]  =  +  57  a^  . 

o  ,  n  n  s'  ,  n 
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the  general  expression  (13)  for  the  variance  of  the  power  is  equal  to 
„  /  N  \  2 


2  A 
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(A6) 


n=  1 


plus  terms  of  the  order  of  1/N  and  higher  compared  to  (A6).  For  convenience 
we  give  (13)  again  here: 
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where,  under  the  assumption  of  real  error-free  fields  used  in  deriving  the  approxi¬ 
mation  (55) 
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Of  the  six  terms  on  the  RHS  of  (13),  the  fourth  and  fifth  are  equal  for  the 
special  case  (Eq.  (A5)]  we  are  considering,  while  the  last  is  of  the  order  of  1/N 
compared  to  the  sum  of  the  fourth  and  fifth  terms  because  of  (Al).  Hence  we  focus 
our  attention  on  the  first  three  terms  of  the  approximation  (13).  The  first  and  third 
terms  combine  to  give 

Now  from  Eq.  (A4) 
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cos  1 2  d  (u  -  us))  =  cos  |  (N  -  1)  m7r] 

so  that,  for  Eq.  (A5)  to  hold,  (N  -  l)m  must  be  an  even  integer.  (Since  m  is  odd 
this  implies  that  for  the  case  we  are  considering  there  are  an  odd  number  of  array- 
elements.  )  But  then 
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As  above  in  proving  the  proportionality  (Al)  we  assume  that  the  amplitudes  of  the 
element  excitation  coefficients  can  be  expressed  as  the  polynomial  in  Eq.  (A2)  and 
hence,  using  Eq.  (A3) 
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(-  1)  n  is  a  polynomial  of  degree  k  in  N  -  1, 
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so  that  for  large  N 
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that  is,  Fq(u)  for  large  N  behaves  as  a  constant  independent  of  N.  Hence,  using 
the  proportionality  (Al) 
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so  that  the  sum  of  the  first  and  third  terms  of  approximation  (13)  is  of  the  order 
1/N  compared  to  the  sum  of  the  fourth  and  fifth  terms. 

As  far  as  the  second  term  of  (13)  is  concerned,  just  as  Eq.  (A9)  follows  from 
Eq.  (A7),  so 
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G  (u)  =  ±  ^  (-  1  )n  "  |w  |3 

n=  1 

and  in  the  same  way  that  we  established  the  approximation  (A  10),  a  similar 

relation  holds  for  G  (u)  with  the  summation  on  k  having  the  upper  limit  31V1  instead 
°  2 

of  M.  It  follows  that  the  second  term  is  of  the  order  1/N  compared  with  the  sum 
of  the  fourth  and  fifth  terms  of  approximation  (13).  Hence,  as  claimed,  for  the 
special  case  we  are  considering,  the  expression  (13)  for  the  variance  of  the  power 
is  given  by  expression  (A6)  plus  terms  of  the  order  1/N  and  higher  compared  to  (A6). 
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